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With the advent of chiral fermion formulations, the simulation of light valence quarks has finally become 
realistic for numerical simulations of lattice QCD. The simulation of light dynamical quarks, however, remains 
one of the major challenges and is still an obstacle to realistic simulations. We attempt to meet this challenge 
using a hybrid combination of Asqtad sea quarks and domain-wall valence quarks. Initial results for the proton 
form factor and the nucleon axial coupling are presented. 



1. INTRODUCTION 

The simulation of light dynamical quarks con- 
stitutes one of the major challenges in contempo- 
rary lattice gauge theory research. The compu- 
tational obstacle has in the past been addressed 
by investigating new algorithms, see, e.g. P] and 
references therein. 

Recently, a new approach has been attempted 
employing a so-called "hybrid" scheme where dif- 
ferent types of sea and valence quarks are used, 
see [2j- It turned out that this calculation is 
extremely successful in predicting spectroscopic 
mass splitting in heavy quark physics correctly. 

In this work we employ another type of hybrid 
calculation involving improved Kogut-Susskind 
sea quarks (Asqtad action [2j) and domain- wall 
valence fermions [!]■ This scheme — although 
breaking unitarity at finite lattice spacing — will 
still have the same continuum limit as a fully 
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dynamical calculation provided this limit exists. 
The quark masses of sea and valence quarks have 
to be properly tuned. The quantities we investi- 
gate are special cases of generalized parton dis- 
tributions (GPDs) 5 . These have both parton 
distributions and form factors as certain limits 
and have already been studied on the lattice us- 
ing conventional schemes [§]■ 

This presentation is organized as follows: After 
discussing the tuning of the parameters of the va- 
lence domain- wall action in Sectional we present 
the main results in Section Finally, we sum- 
marize our findings and give an outlook for our 
ongoing research in Section 0] 

2. CHOICE OF PARAMETERS 

The two essential parameters we have to set are 
the size of the fifth dimension, L5, for the domain- 
wall fermions and the bare quark mass parameter, 
(am)™ F . The tuning of these two parameters is 
discussed in this section. The resulting parame- 
ters and sample sizes are summarized in Tabled 

We used MILC configurations both from the 
NERSC archive and provided directly by the col- 
laboration. We then applied HYP-smearing [7j 
and bisected the lattice in the time direction. 
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Table 1 

Parameters of our runs after the tuning of the DWF valence quark mass. The table shows the lattice 
volume, O, the number of configurations in each sample, the bare Asqtad masses for the sea and 
valence part, the bare domain- wall fermion masses, and the resulting ratio of pseudoscalar to vector 
meson masses. 
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Thus far, we have only considered the first half- 
lattice. We have chosen the gauge field config- 
urations separated by 12 trajectories. In these 
samples we did not find residual autocorrelations. 

2.1. Setting L 5 

The goal in setting L5 is to describe the physics 
adequately at minimal computational cost. For 
finite values of L5 there is a residual explicit chi- 
ral symmetry breaking characterized by a residual 
mass, (am) lcs . We adopt the following definition 
(see jS] for details): 



^Al = 2m q J a q {x) + 2Jl q {x) 
where 

J 5 a 9 (x) «m rra J 5 a (x), 



(1) 



(2) 



which holds up to 0(a 2 ). We require m res to be 
at least one order of magnitude smaller than m q . 

To explore their dependence, we have run sim- 
ulations using two samples of 25 configurations 
with volume Q = 20 3 x 32 from Table □ three 
degenerate dynamical Asqtad quarks with masses 
(a77i)* aqtad ' Boa = 0.050 (denoted as "heavy") and 
two plus one quarks with masses (ara)* aqtad,Bca = 
0.010/0.050 (termed "light"). 

The resulting residual masses obtained from 
Eqs. and J5J are plotted in Figure Q In the 
light quark case, L5 = 16 just fulfills our require- 
ment, while in the heavy quark case L5 = 16 more 
than satisfies it. 

Furthermore, we require the absolute change 
in the masses of the pion and the nucleon under 
changes of L5 to be small; these are shown in Fig- 
ure|3for the pion and nucleon masses. Note, that 



Figure 1. Residual quark mass as a function of 
L5 for the two samples (heavy and light) of 25 
configurations each. 




the errors are not independent between different 
points of L5, since an identical sample of gauge 
field configurations has been used. Therefore, 
we have performed a separate analysis by Jack- 
knifing the differences in the masses compared to 
L5 = 16. The results are shown in Figure 

In the case of heavy quarks, the influence of 
increasing L5 beyond L5 = 16 is negligible in all 
observables. In the case of the light quarks, the 
influence is at most a few percent when going be- 
yond L5 = 16. Hence, we choose L5 = 16 to be 
a good compromise between accuracy and perfor- 
mance. 
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Figure 2. The pion mass with heavy and light 
quarks, m Wl hoavy and m Wl i isbt , and the nucleon 
mass with heavy and light quarks, m Nihoavy and 
m N. n g ht i as a function of L 5 . 
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Figure 3. The same quantities as in Figure|2 but 
with a Jackknife analysis of the mass differences 
vs. L 5 = 16. 
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2.2. Setting the quark mass 

As outlined in the introduction, a necessary 
condition for a hybrid calculation to describe the 
physics of full QCD is that the masses of sea and 
valence quarks are identical. 

Hence, we have chosen to match the pion 
masses of two different calculations: (i) full QCD 
calculations using three and two plus one flavors 
of dynamical Asqtad sea fermions and Asqtad va- 
lence fermions. These results have been taken 
from (ii) Our hybrid calculation with Asqtad 
dynamical sea fermions and valence domain-wall 
fermions with L$ = 16. 

To make the quark masses coincide, we tuned 
the bare mass parameter in the domain-wall va- 
lence fermion action so that the resulting pseu- 
doscalar meson mass coincides with the corre- 
sponding meson mass in the Asqtad case. This 
condition is useful since the pseudoscalar meson 
mass is particularly sensitive to the choice of the 
quark mass. 

The resulting choices for the bare quark masses, 
{am)q, for Asqtad and domain- wall fermions are 
shown in Tabled It is interesting that the dif- 
ference between the bare quark masses for Asq- 
tad and DWF valence quarks is quite substantial. 
Since this difference corresponds to the ratio of 
quark mass renormalization constants, we observe 
that either one or both of these renormalization 
constants deviates substantially from one. 

Finally, we list the resulting nucleon masses — 
which we need for the calculation of the hadronic 
structure — in Table El There is a perceptible 
difference at the heaviest quark mass. The other 
masses are statistically compatible. Furthermore, 
there is no visible finite-size effect at the lightest 
quark mass. Thus, we only observe the expected 
0(a 2 ) effects in the difference between the two 
masses at the heaviest quark mass. 

3. HADRONIC STRUCTURE WITH 
LIGHT QUARKS 

Within the framework discussed in the previous 
sections, we compute different observables of the 
nucleon structure for different quark masses. The 
results are compared to previous results from full 
QCD calculations with heavy quark masses, see 



Table 2 

Nucleon masses as measured with Asqtad and 
domain- wall fermions for the valence quarks. 
SI (ttm)^'" 1 '" (qm)* 3qtad (am)7 F 
20 3 x 32 O050 1.057(5) 1.029(9) 
0.030/0.050 0.930(3) 0.941(10) 
0.010/0.050 0.779(6) 0.756(21) 
28 3 x 32 0.010/0.050 0.763(12) 



References ^U] f° r their calculation. 

3.1. Form factor and transverse size of the 
proton 

First, we consider the electromagnetic form fac- 
tor Fi(—t) of the nucleon. Phcnomcnologically, it 
follows a dipole-shaped behavior. Its derivative at 
the origin specifies the root-mean squared radius, 
r MS , of the transverse charge distribution in the 
infinite-momentum frame and thus characterizes 
the size of of the nucleon. In a hypothetical world 
with heavy pions, the size of the nucleon should 
be smaller due to the absence of a significant pion 
cloud. 

This qualitative behavior is well reproduced by 
our lattice data. Figure 0] shows our first re- 
sults for the proton form factor, Fi(—t). The 
triangles denote the heaviest quark mass in Ta- 
ble 2] the circles the intermediate quark mass of 
(am,)* sqtadsoa = 0.030/0.050 and the boxes the 
lightest quark mass, all at the volume of SI = 
20 3 x 32. The dotted, the dashed, and the dash- 
dotted curves show the best dipole fit to these 
data points. The solid line shows the best dipole 
fit to the experimental data. 

Figure [S] shows the transverse mean square ra- 
dius, r MS . The boxed data points represent the 
old SESAM data [Ej, obtained in a full QCD 
calculation with Wilson quarks. The triangu- 
lar data points show the hybrid calculation with 
the data points from Table with the volume 
SI = 20 3 x 32. Finally, the circle represents the 
volume SI = 28 3 x 32 in Table [IJ As the quark 
mass decreases, the transverse size of the nu- 
cleon increases and approaches the experimental 
value. Note that we find no evidence of substan- 
tial finite-size corrections since there is no sta- 
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Figure 4. The proton form factor, Fi(— t), as a 
function of the virtuality, t = q 2 . 
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Figure 5. Transverse root-mean squared radius 
r MS of the nucleon as a function of the pion mass 
squared. 
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tistically significant difference between the data 
points at different volumes. 

3.2. Axial charge 

The axial charge, gA, is given by the forward 
value of the generalized parton distribution ^4^ d 
(see Reference [B] for details). This quantity 
has been shown to exhibit substantial finite-size 
dependence an d therefore provides a good 
benchmark of the importance of the physical box 
size as we approach the chiral limit. 

Figure [fj] shows a central result of this work, 



the nucleon axial charge as a function of the pion 
mass squared. It has been renormalized using the 
five-dimensional conserved current. The symbols 
are the same as in Figure [S] 



Figure 6. Nucleon axial charge, gA, as a function 
of the pion mass squared. 
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We observe that the SESAM results in a 
(1.6 fm) 3 box first increase with decreasing pion 
mass and then fall off when the pion cloud be- 
comes too large to fit in the box. In the larger 
(2.6 fm) 3 MILC box (on a 20 3 lattice), the value 
at nift ~ 700 MeV increases back to its physi- 
cal value again and then falls off at ~ 300 
MeV when the light pion no longer fits in the 
box. Increasing the box to (3.5 fm) 3 (28 3 lattice) 
allows gA to increase to a physical value consis- 
tent with experiment. Thus, the locus of points in 
the largest box at each mass extrapolate smoothly 
to the experimental result. 

To summarize, we give the resulting values of 
r MS and gA in Table |21 

4. SUMMARY AND OUTLOOK 

In this talk we have used a tuned hybrid cal- 
culation of Asqtad sea quarks and domain-wall 
valence quarks. We have determined an optimal 
value for L5 in the domain-wall action and tuned 
the bare quark masses such that the pion masses 
of an Asqtad- valence and the domain-wall valence 
calculation agree. If the continuum limit exists, 
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Table 3 

Summary of the hadron structure results. The 
table shows the volume, the Asqtad sea quark 
mass from Table ^ an d the resulting values of 
r MS and g A . 
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this calculation provides a valid scheme to com- 
pute all hadronic observables. 

We have computed the F\{— t) form factor of 
the proton and obtained the r MS radius from it. 
The results are in qualitative agreement with the 
picture of the nucleon becoming larger as the pion 
cloud becomes more prominent with decreasing 
quark mass. 

We have also applied our scheme to the case of 
the nucleon axial coupling, g A - Whereas calcu- 
lations in a fixed volume underestimate g A when 
the pion becomes too light to fit in the box, our 
sequence of calculations in three volumes produce 
a locus of points that smoothly extrapolate to the 
experimental result. 

Encouraged by the successes reported here, we 
are continuing the calculation of hadronic observ- 
ables using improved staggered sea quarks and 
domain- wall valence quarks. 
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